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ABSTRACT: 


With the assistance of two extra groups, ie., an extra hidden gauge group SU(2)p 
and a global U(1) group, we propose a two component dark matter (DM) model. After 
the symmetry SU(2)p x U(1) being broken, we obtain both the vector and scalar DM 
candidates. The two DM candidates communicate with the standard model (SM) via three 
Higgs as multi-Higgs portals. The three Higgs are mixing states of the SM Higgs, the 
Higgs of the hidden sector and real part of a supplement complex scalar singlet. We study 
relic density and direct detection of DM in three scenarios. The resonance behaviors and 
interplay between the two component DM candidates are represented through investigating 
of the relic density in the parameter spaces of the two DMs masses. The electroweak 
precision parameters constrains the two Higgs portals couplings (Am and 62). The relevant 
vacuum stability and naturalness problem in the parameter space of Am and ô> are studied 
as well. The model could alleviate these two problems in some parameter spaces under the 
constraints of electroweak precision observables and Higgs indirect search. 
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1 Introduction 


The discovery of Higgs particle with mass around 126 GeV at Large hadron collider (LHC) [1, 
2], roughly consistent with the SM predictions, seems complete the SM. However, the nat- 
uralness problem is still unsolved. To alleviate the problem, as well known and extensively 
studied, more bosonic fields are needed [3-7]. Recently, it has been noticed that the Higgs 
field-strength renormalization could be enhanced by these new boson fields [8-10]. And 


that after the Higgs field normalization, the Higgs coupling's modification could be de- 
tected indirectly [9, 11]. Thus, it shed light on testing the mechanisms which alleviates the 
naturalness problem. Secondly, the SM has to be extended in order to accommodate the 
cold dark matter (CDM) and baryon number asymmetry of the universe (BAU). In order 
that the BAU not being washed out after its generation, the strong first order electroweak 
phase transition (SFOEWPT) is necessary, which is out of the capacity of the SM and new 
scalar fields are needed [12-14]'. At the same time, the way to test the mechanism allevi- 
ating the naturalness problem could also been used to test SFOEWPT, which is just one of 
the preliminary aims of the International Linear Collider (ILC) [8]. Thirdly, vacuum insta- 
bility could be rescued by introducing vector and scalar fields [3, 17-22], and thus making 
the inflation with not very small top quark mass possible to coincide with the BICEP and 
Plank data [23, 24]. 


All above arguments require us to extend the SM with new bosonic freedoms, which 
we choose vector and scalar fields. There are many studies on the multi-component DM 
scenarios, see |3, 25-35], as well as the dynamical DM scenario [36] whose phenomenological 
consequences are often quite distinct and can be applied to a much broader variety of multi- 
component DM scenarios [37]. The primary purpose of this work is to explain the CDM relic 
density, where the two different DM components interact with each other besides with the 
three Higgses, thus affects the evolution of the DM components number densities through 
the coupled Boltzman equations [3, 38]. 


For the method to introduce vector dark matter through the effective Lagarangian 
method yielding strongly constraints on the parameters space from the unitarity con- 
straints [3, 39, 40], we consider the scenario in which the vector dark matter fields (V^) 
respects an extra non-Abelian gauge symmetry SU(2)p. After the SU(2) p being broken via 
the the complex doublet(¢), one SO(3) symmetry is induced which V^ respects to, thus 
making V^ stable. One more complex scalar singlet, S, is supplemented to the model. 
After the global U(1) symmetry, i.e., S — e'^S respected by V (H,¢,S), being broken 
spontaneously and softly, the real part of S (S) gets a vacuum expectation value(VEV), 
and the imaginary part of S (A) respects the reduced Z2 symmetry, which makes the other 
DM candidate. After the breaking of the SU (2) p x U(1), the Higgs field of the SM mixed 
with two scalar fields, the real parts of $ (7) and S. Since the two component DMs inter- 
act with the SM particles and each other through the three Higgs fields, one could expect 
three resonance enhancement effects. And three significant enlargement of the magnitude 
of the annihilation cross sections, of the DM to DM and the DM to SM particles, appears. 
'The magnitudes of the relic density around the three resonances decrease, since which is 
robustly inversely proportional to the annihilation cross sections of the DM to SM particles. 
The magnitudes of annihilation cross sections of scalar DM fields (vector DM fields) to the 
SM particles, around the three Higgs field masses, are enlarged due to the t,u and seagull 


! In fact, additional heavy fermions could also make the SPOEWPT feasible, see [15, 16]. 
?The hidden gauge theory has also been used to study self-interacting dark matter [42, 43], inflation [44], 
and the model with a kinetic mixing portal between the gauge boson DM and SM [41]. 


channels of AA(VV) — h;h;®. The spin-independent (SI) DM-nucleon scattering cross 
sections are determined by t-channel interactions between the nucleon and two component 
DMs through the exchangeing of three Higgs fields. The model could be distinguished 
from the model with no interaction between the two DMs, since the opening of the channel 
AA — VV(VV — AA) does affect the evolution of DM number density, and thus the 
magnitude of relic density. This kind of effects are explored and illustrated in DM relic 
density analysis section of this work. The Higgs indirect search at the LHC requires small 
mixing between h and y (S). The electroweak precision observable experiment imposes 
stringent constraints on our parameter spaces. The vacuum stability of the model could 
be improved under the above two considerations. We explored the way to alleviate the 
naturalness problem and its’ indirect search as well. 


The paper is organized as follows. We construct the model and explore relic density 
and direct detection of the DMs. After which, we investigate the Higgs indirect search and 
electroweak precision constraints on the model. And then the stability of the model are 
studied. We also consider the possibility to alleviate naturalness problem and the way to 
trace the footprint of which. At last, we conclude this work with discussions and conclusion. 


2 The Model 


To construct the model which includes the stable vector and scalar fields, we introduce 
these two fields as follows. The vector field V is introduced as the gauged field of the 
SU(2)p symmetry, which couples to the SM through a doublet, ¢, which is a singlet of the 
SM but charged under the SU(2)p. Here, one should note that the mixing between V'^ 
and the SM gauge bosons through kinetic mixing is absent for the non-abelian character 
of the SU(2)p. We supplement one more complex singlet S. And one more global U(1) 
symmetry is required besides the SU(2)p and the SM group. The complex singlet, which 
interacts with the SM and the SU(2) p through the Higgs portal and ¢ portal terms in the 
scalar potential, transforms trivially under the SM and SU(2)p gauge group. 


The Lagrangian of our model is 
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with 


Vis = Vno + Vus + Vos, (2.2) 


?Hereafter, the Lorenz index yp of the vector field V^ will be hided for simplicity in some cases. 


and 


Vio = Amb OHH, (2.3) 

Vins = 2H HISP, (2.4) 
ô 

Vos = 59 lS? , (2.5) 


where Do = Oh — ir - V'^. The hidden gauge coupling gg < 47, required by the 
unitarity bound, need to be hold for any thermal particle whose relic density arises from 
the freeze-out of its annihilation [39, 45, 46]. In particular, bı and a, terms break the global 
U (1) symmetry explicitly. In the SM Lagrangian, the Higgs potential notations are defined 
as: £?M > —(m?/2)H'H — A(HTH)? with H = (0,v + h)/V2, where v is the VEV of the 
Higgs field. 


After the SU (2)p being spontaneously broken, both ¢ and the singlet S get the VEVs: 
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So the Eq. (2.1) in the unitary gauge recast as, 
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cFVhawy,S,A , (2.8) 
here V^ = UV'^U-1— i[o,U]U with U = exp(—it -€/vg). The vector DM mass is given 
by my = ggvq/2, and the tree-level potential Vo (h, m, S, A) recast as, 
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Vo (h, n, S, A) = UI 


+2 ((S + v.)* +A’) . (2.9) 
Here, we would like to mention that, 7/ lives in the fundamental representation of SU (2) p, 
and displays a custodial symmetry SO(3) in the Vr 9,3 component space, which makes three 
V” components degenerate in mass and thus stable [47]. The explicit Zo-breaking term 
is proportional to a4 and being introduced here to avoid the cosmological domain wall 
problem [48-50]. After choosing ¢a, = b, = 7, the potential retains a Zo symmetry for 
Im(S), thereby ensuring the stability of the particle A [51, 52]. 


Requiring that the potential in Eq. (2.9) has a minimum at (H) = h/V/2 = 0 and 


(S) = S +iA=0 +i- 0, the following minimization conditions are obtained: 
OVo OVo OVo OVo 
E - l = = 2.1 
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where all derivatives are evaluated at (h, n, S, A) = (0,0,0,0). These minimization con- 
ditions allow the Higgs VEV v and the singlet VEV v, to replace m? and by according 


to 
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Thus, at the minima, the mass matrix is obtained as 
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And in the basis of (h, 7’, S), we have 
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To work in the mass eigenstates, i.e., h1,2,3, we diagonalise the mass matrix Eq. (2.13) 
through 


RM? RT = Mi, (2.14) 
with matrix R being given by 
C103 C381 $3 
R- C281 — €18283 C1C2 — 815253 C359 ; (2.15) 


$182 — €1€253 C152 — €25153 C2C3 


c, s and their subscripts 1, 2, 3 represent cos, sin, 015, 053, and 043 individually. 


3 Dark matter analysis 


In our model, we have two component DMs. Thus, to present the novelty of the model 
properly, we use the coupled boltzman equations explored in [3], with annihilation channels 
being depicted in Fig. 1, and details of annihilation cross sections are listed in section A. 
The relic density of each component DM could be obtained through Q4 yh? = 2.755 x 


108 MAY. Ya v(To) [53, 54| after we calculated YA v(To) numerically, and the total relic 
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Figure 1. Feynman diagrams of annihilation channels. 


density is the sum of the two component DMs, Qh? = Q4h? + Qyh?. As for the numerical 


calculations of YA y (To), we refer to Eqs. (B.1), (B.2), (B.3), and (B.4). 


In order to investigate how does the interplay between the two DM components affect 


the evolution of the DM abundance and the dependence of the DM relic density on each 


parameters, we scan the parameter spaces according to the following three groups: 


1. Scan in the M4 — My plane for the case of M4 > My, with the other parameters 


being given in Table 1. 


2. Scan in the Ma — My plane for the case of M4 < My, with the other parameters 


being given in Table 1. 


3. Scan in the Àm — 62 plane, with DM masses being fixed as: My=120 GeV, M4—150 


GeV, and the other parameters are given in Table 2. 


U Ug | Us À Ag | do 1 02 Am Jo 
246 | 738 | 123 | 0.515 | 0.2 | 0.2 | 0.02 | 0.04 | 0.01 | 0.04 


Table 1. The input parameters for the analyses of the cases of m4 > my and mA < my. 


U Ud Us AÀ Àg d» ôi 
246 | 738 | 123 | 0.515 | 0.2 | 0.21 | 0.02 


Table 2. The VEVs of three scalar fields and other couplings. 


3.1 Relic density analysis 


Set free parameters as in Table 1, the three Higgs masses are obtained as mj, = 124.8, mp, = 


233.8, and mj, = 71.1 GeV in order. 


1. The first group with mA > my 
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Figure 2. Left: Plots for the cross sections of channels AA —^ X X and AA — VV in units of 


0.3894 x 10? pb; Right: The corresponding relic density Q4A?. 
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Figure 3. Left: Plot of the cross sections of channel VV — X X in units of 0.3894 x 10? pb; Right: 
Corresponding relic density Nyh? for m4 > my. 


The annihilation channel AA — VV is opened up and shown by the green line on 
the left panel of Fig. 2. The bandwidth of the green line comes from the change of the 
vector DM mass my, and the bandwidth becomes wider and wider with the increasing 
of scalar DM mass ma, since the bigger the scalar DM mass mA the larger viable 
range will be left for my. Fig. 2 depicts that the channel AA > X X is the dominate 
one when mA >70 GeV, and the variation tendency of the relic density (4A? is plotted 
in the right panel. Three peaks of the annihilation cross section of AA — VV are 
presented, at m4 = Mh, /2, Mha /2, Mha /2, those are the three resonances at these 
mass values, as could be seen from Eq. (A.6) as well. At the same time, we find 
three other peaks of the annihilation cross section of AA + XX, besides the same 
three peaks as that of AA > VV, at m4 = mp, (the three masses of three Higgs), 
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because the new partial t, u and seagull channels AA > hih; through exchanging the 
scalar DM particle themselves are opened up around these mass values. These peaks 
represent themselves on the right panel of Fig. 2 by the transition of the color along 
axis of m4. When there is a peak bigger than about 1 x 10-? for the cross section 
(on the left panel), a decrease behavior of the corresponding relic density appears on 
the right panel. 


Last but not the least, we would like to emphase that the opening of the channel 
AA — VV does cause the decrease of the magnitude of Q4h?. Because of Q4h? ~ 
1/(ov) 4a-+xx, when the channel AA — VV has been shut down, one may expect 
(cv)AA S xx ~ pb and Q4Rh? ~ 0.1. While, the right panel of Fig. 2 depicts more 
decrease of the magnitude of Q4h? around m4 = Mhi 5 /2, which demonstrates the 
effects of the annihilation process AA — VV. 


For the case of m4 > my, as shown in Fig. 3, the channel of VV — AA is closed 
when we analyzing Qyh? . The cross section of the channel AA + VV together with 
that of VV — X X affect the corresponding relic density Qyh?, and the second one 
dominates the value of Qyh?. Only the slim peak caused by the resonance at about 
my = mp, /2 brings a sizable decreasing of Qyh?. The viable region of my matching 
the experimental value Qh? = 0.1189 [55|(Qyvh?) is very small as shown in the right 
panel. 


2. The second group with my > mA 
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Figure 4. my > ma. Left: Plots for annihilation cross section of channel AA > X X in units of 
0.3894 x 10? pb; Right: Relic density Q.4h?. 


In this case, the channel AA > VV is forbidden due to my > ma, and the behavior 
of the cross section of the channel AA > XX is given in the left panel of Fig. 4. 
With the increasing of m4, the first peak caused by the resonance effect does not 
leave any trails, as could be seen in the magnitude of the corresponding relic density 
(see the right panel of Fig. 4), because the magnitude of the cross section is too 
small. The second peak is caused by the resonance effect at m4 = mp,/2, and the 
tinny increasing of cross section induces a small decreasing of the corresponding relic 
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Figure 5. my > mA. Left: Plots for the cross sections of annihilation channels VV —^ X X and 
VV — AA in units of 0.3894 x 10? pb; Right: The dominating relic density Qy ^?. 


density. The third peak demonstrates the opening up of the partial t and u channels 
for AA — haha and also brings a very big decrease of the relic density at about 
mA > Mp, as expected. And the fourth peak, which is caused by the resonance at 
MA = Mh, /2, brings a sizable decrease of the magnitude of the relic density also. We 
also notice the small waves at about my = 160 GeV in the right panel of Fig. 4. The 
phenomenon could be explained by Eq. (B.4), which illustrates that the behavior of 
the cross section (ov)vy_5 AA affects the calculation of Q Ah? and the variation of the 
magnitude of the cross section (ov)vy., AA generates these small waves. The channel 
VV — AA also causes a bigger decrease of the magnitude of Qy A? in comparison 
with the scenario in which the channel VV — AA has been shut down, as shown in 
the right panel of Fig. 5, especially around my ~ mp, /2(mj,). 


Fig. 5 is plotted to show the magnitude of Qy A? contributing from both channels 
VV > XX and VV — AA. The cross sections are small for most regions, and a big 
decrease of Qyh? exists in the two regions my ~ Mh, /2 and my > Mp. For the first 
decrease, both channels induce comparable effects. And for the second decrease, it is 
the channel VV — XX that dominates the magnitude of Qyh?. Especially when 
my > Mh, the partial t, u and seagull channels VV — hoho through exchanging the 
vector DM particle are opened up. These effects give rise to too big (ov)vv., xx, 
and then the door to generate relic density Qy A? is almost closed. 


At last, the two panels in the Fig. 6 presents the total relic density for m4 > my 
and my ma, respectively. Under the set of parameters as in the Table 1, the large 
enough magnitude of the relic density could be obtained easily. 


3. The third group: varying Àm and 6 


FIG. 7 depicts the relationships among 62, Àm and the cross sections. FIG. 8 illustrates 
that the magnitude of Q4AA? (Qv?) reaches its critical value around 62 ~ 0.09 (Am ~ 
0.01). And that the DM relic density for each component almost depends on their 
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Figure 6. Plots of the total relic density 04h? + Qy h? for the case of m4 > my and my > ma. 
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Figure 7. Plots of ô and Am vs. cross sections. 


own coupling parameters, this is caused by the fact that all related annihilation cross 
sections are proportional to d2 or Àm. 


3.2 Direct detection 


In our model, the two component DMs, i.e., the scalar A and vector V, interact with the 
SM particles through the exchange of three Higgs bosons. Thus, the DM-nucleon scattering 
cross section is spin-independent (SI). For each single component, the SI DM-nucleon cross 
sections are calculated to be 
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with Rij coming from Eq. (2.15), my and fy = »; fL -3x xu are the nucleon mass 
and effective Higgs-nucleon coupling respectively, where fy is the summation of light quark 
(fr) and heavy quark (fg) contributions, and we take the value fw = 0.326 [56] in our 
numerical analysis. Since the current experiments assume that the local DM density is 
provided by one single DM specie, the situation that both A, V components contribute 
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Figure 8. Plots of 05 and Am vs. relic density. 


to the local DM density making us unable to use the current experimental results directly. 
Assuming the contribution of each DM component to the local density is the same as their 
contribution to the relic density, the SI scattering cross section should be rescaled by a 
factor Q4 y h?/Q py h?. Thus, the corresponding upper limit on the SI DM-nucleon cross 
section of each single component is [38, 52]: 


F5 (omh /Qsh?)e ST (MA), (3.3) 
oS € (pu? /Qyh? oS? (Mv). 


Similar to the above Section, here we do analyses using the LUX experiment results [57] 
via three groups as well. 


1. mA» my 
The left and right panels of Fig. 9 depict that: with the value of my being unbounded, 
ma should be bigger than 30 GeV. 

2. my > mA 


The left and right panels of Fig. 10 depict that mA,y should be no smaller than 
50 GeV. 
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Figure 9. For the case of m4 > my, parameter regions being left considering the direct detection 


constraint of LUX. 
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Figure 10. For the case of my > ma, parameter regions being left considering the direct detection 


constraint of LUX. 
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Figure 11. 62 and Am. Parameter regions being left considering the direct detection constraints 


on the model coefficients. 


3. parameter space of ó9 and A; 
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mA 


mA 


Fig. 11 depicts that all the parameter spaces satisfy the experimental constraints. 


4 Higgs indirect search and electroweak precision constraints 


In this section, we study the effects of the mixing among the three fields h, r/, S using the 
Higgs indirect search and constraints from electroweak precision observables. 


4.1 Higgs indirect search 


'The current status of the LHC measurements of the Higgs couplings constrains the matrix 
element R11, which describes the discrepancy between hı and the SM Higgs. The couplings 
of h4 to all SM particles are the rescaled values of the SM couplings, taking the form of 


gn xx = Rugs xx - (4.1) 


For the first two group parameter spaces in the last Section, additional decay channel 
hı + AA(VV) exists, which almost does not change the total width, since we have the 
square of R15(R33) to suppress the magnitude of P(h; + AA(VV)). The same logic applies 
to the analysis with the third parameter group used in the last Section +. Thus signal rates 
Lxx associated with Higgs measurements are functions of R11 


o: BR 


_ E 
LXX = SW. BRIM ^ Rh; (4.2) 


with c, BR (that with a a superscript S M) being the production cross section and branch- 
ing ratios of hı (the SM Higgs). Therefore, to what extent the model differs from the SM 
Higgs measurements is determined by the value of R41. The value of R11 for the parameter 
setup given in Table. 1 is 0.9973, and we refer to Fig. 12 for the R vales? corresponding to 
the benchmark scenario given in Table. 2. 


4.2 Electroweak precision observables constraints 


The presence of the scalar hg with mass mj, < 114 GeV is subject to tight bounds from 
the LEP [58], while these bounds are highly released in all our parameter choices, since 
the mixing between hı and hg is very small, about 0.0573 for benchmark scenario given in 
Table. 1 and could be seen in the Fig. 12 for the benchmark scenario given in Table. 2. 


As for electroweak precision observables, the oblique parameters S and T' which pa- 


rameterize potential new physics contributions to electroweak radiative corrections, are 


“Because the mass of ha(h3) is above the threshold, additional channel hi — hzh2(h3h3), which changes 
the total width, does not open. For the study of one Higgs decays to two other Higgs bosons, we refer 
to [59]. 

?Here, we use R(i, j) as the labels in the two plots, which means the same as that of Ri; in other place 
of this work. 
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Figure 12. R values with respect to Àm and 62. 
computed following Refs. [47, 60]. In our model, they are 


1 
S= me [5 [log Ti. + G(m;, , mZ) m G(m;, , mZ)] 


+Ri2 [log Rash T G(m;,,, mz) E G(m;,,. mz)| 


TRI, [log Rhah + G(m;,, , m7) — G(m;,,. mz)| H (4.3) 


T= 3 R2 | 1 ( log Rzn, log Rzn ) ( log Rwn, log Rwn ) 
167 sin? Ow 1l | cos? Ow X1— Rzn, 1— Rza 1— Rwn, 1— Rwn 
1 log RZhə log Rz log Rwho log Rwn 
cos? Ow \1— RzZho 1— Rza 1— Rwn, 1— Rwn 


+R? 1 log RZhs log Rz log Hwna log Rwn (4 4) 
P cos? Ow 1— RZhs 1— Rz 1— Hwna 1— Rwn i l 


+R? | 


where, Rap, G(m%,m%) and f(RAp) are given by 
mA 


RAB = —2. 5 (4.5) 
mB 


79 
G(mA,m2j)) = —— + 9Rap — 2RÀg + (12 — ARAg + Rig) f(RaB) 


3 
R 1 
+(—10 + en — 6Rå4p + eg s i ) log RAB, (4.6) 
Rap-l 
Ran hip c d) ZT log RaB—2-y dune Rap > 4, 


24/Rag(4-— Rap) arctan 4/ fae Rap « 4. 


The parameters set is given in Table 2 and the magnitude of gy is determined by my 
and vg. The constraints on two parameters 62 and Àm from S and T parameters are shown 


— 14- 


Figure 13. Constraints on à5, Am from limits of S and T. 


in Fig. 13. The S and T constraints require that 62 and Am should be smaller than about 
0.45 and 0.06 respectively. S and T' are very sensitive to the mixing effects. The first 
two parameter groups (as given in Section 3.1) does not change the mixing angles among 
h, ņ', S, and S and T give rise to the null limits since the mixing effects among h and 7, S 
in that cases are very small. 


5 Vacuum stability 


The global minimum of the tree-level potential of our model requires, 
A> 0,Ag > 0, dp >0, Adz > 402 , Agdz > 02 , Aàg > 16A, . (5.1) 


From one-loop renormalization group equations (RGEs) of the Higgs quartic couplings in 
Appendix C, i.e., Eq. (C.1), we find that the Higgs portals couplings Am and 42 are all 
get involved and give rise to positive contributions to £4. Thus, we could expect the 
vacuum stability problem being solved or alleviated in some parameter spaces. Adopting 
the central values of top quark mass, the Higgs mass, and the strong coupling [55] as the 
low energy boundary conditions, we find that with the increasing energy scale, the Higgs 
quartic coupling running to a negative value around the scale 10? GeV, and then grows to 
positive values latter |61, 62]. 


Based on arguments on vacuum stability given in Ref. [62], one needs have positive 
value of the Higgs quartic coupling to ensure absolute stability. It is obvious that to obtain 
absolute stability, we should elevate the curves of the Higgs quartic coupling in the plot of 
A — n |61, 62]. From the 8 functions given in Eq. (C.1), we find that with the increasing 
of Àm and à» one could have the increasing of the value of the Higgs quartic coupling A. 
To verify this, we explore parameter spaces that survive under S and T' limits as shown in 
Fig. 14, and we find that Am and 62 are required to be bigger than around 0.35 in order to 
evade the stability problem. 


To figure out to what extent we need to lift the value of the Higgs quartic coupling in 
the third parameters setup (as given in Table 2), we plot contours of A(u = 10? (1019) 
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Figure 14. The Higgs quartic coupling A as a function of Am and d2 at 10? GeV (left panel) and 
1015 GeV (right panel). With parameters set as: A = 0.129, Ay = 0.2, dp = 0.2, 04 = 0.2, dp = 0.2, 
às 203. 


GeV) with respect to Am and 62, as shown in Fig. 15, and we find that the vacuum is not 
bound from bellow there. The plot of A (at the scale of 10? GeV) as a function of Am and 
62, for the scenario with two S being supplemented (both two S have the same interaction 
with the SM and the SU(2)p group), shows that with 62 > 0.25 one can indeed obtain 
the potential bounded from below, as shown in Fig. 16. The two figures, i.e., Fig. 15 and 
Fig. 16, illustrate that to achieve the vacuum stability, no smaller than two S is expected 
in the third parameters setup. 
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Figure 15. The Higgs quartic coupling A as a function of Àm and d2 at 10? GeV (left panel) 
and 10!5 GeV (right panel). With parameters set as: A = 0.515, Ay = 0.2, d2 = 0.2, 6; = 0.02, 
52 = 0.04, Àm = 0.01 . 


6 Footprint of the naturalness problem 


In this section, we would like to present the relevant naturalness problem in our model, and 
discuss the indirect search of the scenario which could improve the naturalness problem. 
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Figure 16. A(10? GeV) as a function of 62 and Am. n = 2 is to indicate we have supplemented two 
S. 


6.1 The naturalness problem 


In the SM, considering the gauge invariant property of the two point Higgs Green func- 
tion [63], the naturalness problem could be defined as 


A 
(mg) = mh + 7 —;VCsw ; (6.1) 
T 
where 
92,92 
VCsm = 12A — 12g; + 292 A 391 ; (6.2) 
and mp (m9) is the renormalized (bare) Higgs mass, A indicates the cut off scale where the 


new fields are required to cancel the quadratic divergences of the Higgs mass square. 


Although our model is renormalizable, we assume it as a low energy effective theory of 
some more fundamental theory, which is UV completed at the scale A. Thus, Eq. (6.2) 
becomes 

042 2 A? 
(mg) — my d amp OVA , (6.3) 


with 


PES 
VOvA ~ 12A + 592 d gig? 


2 
5 Aj 392m | ,3,. ÖzÀm 
+ ^m LEUTE n( 542+ 7 
dd» | 530 
24> A 


in which new fields get involved. Here, we suppose that not only one S embracing the same 
property as explored in the model (thus n is introduced to denote the number of S), and 


EX m 


the symbol z is to indicate that Eq. (6.3) could be hold in the case where S, n’, and H 
have negligible mixing. 


In addition, we would like to mention that, the VEV and Higgs mass in the SM are gauge 
invariant considering the renormalization of the mass term in the SM Higgs potential [64]. 
And [64] use the different tadpole renormalization method in comparison with that of [63] 
and leaves no trails of tadpole contributions in VCsm and VCy4, which cast the new 


formula of 
I 9 2 3 2 
VCsm = —(6A— 6gt + 793 + 791) ; 
4 4 
9 3 1 n 
VCya © —(6A+ 792 + 391 — 69) — 55m — 392 - (6.5) 


To soften the naturalness problem |4, 5], the simplest solution is to suppose VCy 4 = 0 
(or VCt,, = 0) so that the modified Veltman condition [65] is realized. For A = 10 TeV, 
n = 1(6) is needed to obtain VCy 4 = 0 (VCy 4 = 0). 


6.2 Indirect search for the scenario which alleviates the naturalness problem 


Considering My 5,4 >> v could been satisfied in some parameter spaces of the model 6. one 
may integrate out the n’, S, A and express their effects in terms of an effective Lagrangian 
below the scale Min(m,,mg,m,), which involves only the SM fields with appropriate 
higher-dimensional operators. At one-loop level, integrating out the 7’, S, A leads to shifts 
in the wave-function renormalization and the potential of the Higgs doublet H, as well 
as operators with dimension six and higher. The dimension-six operators in the effective 


Lagrangian cast the form of 


Cie + (cH a, Hy (19 pon ap (6.6) 
eff m SM m2, l m2 m, 2 H . . 
n 
Matching to the full theory at the scale my, 4, we have 
/ Thy! 
di = "lm l?/(96x?), 


TL 
cj = gol / (9677), 


cA = a s / (967). (6.7) 


Below the scale of electroweak symmetry breaking, Eq. (6.6) leads to a shift in the wave- 
function renormalization of the physical scalar h [9], with 6Zp, = 2v?(c7/ m, + c$/ m2 + 
ch /m?,). After canonically normalizing h, i.e., h > (1-6Z),/2)h, its coupling to vectors and 
fermions are altered, which may lead to a measurable correction to, e.g., the hZ associated 
production cross-section 


do zy, = -2v (cd [m2 + cf, /m + cf / mA) , (6.8) 


°For simplicity, we consider the case in which the mixings among the three fields n’, S, A could be 
neglected. 
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where we have defined daz; as the fractional change in the associated production cross 
section relative to the SM prediction, which by design vanishes for the case of the SM. 7 
The sizable ny (n), being required to relax the naturalness problem ? (see Eq. (6.3)), may 
causes the correspondingly observable effect in the precision measurement of ozp, which is 
what the future lepton colliders supposed to detect. At last, if one use our model to analyze 
electroweak phase transition, the measurement of associated production cross section of Zh 


(ozn) might imposes very strong constraint on the multi-Higgs portals couplings A», and 


5. 


7 Discussions and conclusion 


In our model, the vector field (V) and the imaginary part of complex singlet (A) are 
stable due to the reduced custodial symmetry SO(3) and the residual Z2 symmetry. The 
interactions between the two DMs and resonant effects mediated by three Higgs portals 
are demonstrated through the study of the relic density behavior with respect to m, and 
my. The effects of the multi-Higgs portals couplings on the generation of DM relic density 
are shown in the parameter space of Am and 62. The parameter spaces of Am and ô> are 
totally free under the LUX experiment results limits. Constraints (from LUX experimental 
results) on parameter spaces of m4 — my give rise to m4 > 30 GeV (ma(my) > 50 
GeV) for the case of m4 > my (ma < my), which allows more relaxed parameter spaces 
in compare with the simple models. As for the third parameter group: the electroweak 
precision observables require Àm < 0.06 and 62 < 0.45, the vacuum stability at high scale 
requires at least two complex singlet, which is supposed to be consistent with the naturalness 
problem argument. Other benchmark scenario with one S, which survives under indirect 
Higgs search and electroweak precision test limits, and could improve the stability problem 
is also explored . 


At last, we would like to mention that, through tuning of the parameter ô> associate 


with tuning of v, in the first setup of parameter spaces ? 


, one could use the model to explain 
Galactic Center Excess observed by the Fermi Telescope, with (ov) z 107% em3/s being 
provided by scalar DM pairs annihilating to bb and the relic density supplied by the vector 
DM”. Details of which, and the realization of SFOEWPT and inflation in the model are 


left for further studies. 


"Here, we would like to mention that, which one of Ch. CH, cá in the our model takes part in ozn 
depends on parameter choices, i.e, the derivations of et, ch, cà require Tr s, A È v. 

8 The sizable n is also necessary from the viewpoint of the vacuum stability, see Section 5. 

?In order to maintain small mixing between h and S which helps to escape constraints of LHC and 
electroweak observable experiments. 

10For a viable annihilating multi-component dark matter model consisting of two real gauge singlet 
scalars, one can explain the low energy (1— 3 GeV) gamma ray excess from both Galactic Centre and Fermi 
Bubble [67]. 
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A Annihilation Cross Sections 


When the Higgs mass mp, is larger than twice of the SM particle masses, the corresponding 
visible decay channels widths are 


N.G pm? 


4m? i 

Taff = — ma Ral (1 - = ; (A.1) 

Av2n Mp; 

4m? 4m? 12m4 

Th, R Wi Ws Wl, A.2 
hioww = = —— mi, | al? m? | mi. mj, ( ) 

Gr 3 9 4m?, 4m? | 12m? 

= —— m? |R; 1 1 | A.3 

hy ZZ T il m2 zx mi (A.3) 


and when the Higgs mass mp, is larger than the twice of the DM mass, the invisible decay 
channels widths are 


Leno 2 2 4 
"i v Th. Amy, Amy, | 12my, 
T =|R : 1 ] A.4 
mad dm 285 1281 mt my. ms, k Mp, ( ) 
bo, 8 dw v Am? 
U U U m 
Poem Hn 5 t Ra, Ra su A, (A.5) 
2 v| 327mp, m, 


The total decay widthes of the Higgs with mass mp, are the summations of visible and 
invisible decay widths 


'The annihilation cross sections for DM corresponding to Fig. 1 are given as follows 


e Top-left 
4,2 2 
1 99U$ Mp 
< ov >} "x : 
VK hgh 327 mi mi 
Y Ro; (Ra fa; + Fijo Raj Amv + Riz Roj% vs) 
x 2 = 9 P , 
: Amy, — my. + imp, Tn, 
2 2 \ 3/2 2 2 
m m (g5vo/v) Roi Ri 
ES f f $ 
< ov "wan = «Ys i-i) 2 
f i 


2 . 
ms + imp, Uh, 
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rds mi m3 1 m? 3m4 Y (g5vo/ v) Ra Ri 
ov = 
VV>Z2 — 1287 mi mi mi : 4m, — m, imp Up, 
2 R, 2 
< ov >yvoww = my 1 miy 1 K UN B i 
647 mi, mi, mi, : Am?, — my. +imn Tn, 
1 1 m? 
< ov >? PA 1 
TU ZAA>hjh; 7 39r m?, m2, 
is Ry (Ra fij + Fijo Raj Amv + Riz R3; Bus) 
x |v : 
2 2 4m, = m, + imp, Un, 
ôi Roi (Ra Rijrv + Riz Roj Amv + Ri3R3j Bus) 
+ =v : 
2 $ 2 4m, L m; + imp Up, 
do Fi (a fij + Fi Fio; Am Ue T Riz Raj Bus) Á 
+ =y. , 
97 > 4m2, — m, impp, 
» mj 1 m; is 5 (05/2) Rui Ri 
< 0U >AA =n = =n - 
>S ° 7 167 m2, : 4m, — m; impp, 
2 
Y 3 (01v9/2v) Roi Ri y (dov, /2v) Rai Ri 
z 4m2, — m; c impp, : 4m, — m; - impp, 
m2, m? m? 3m4 (02/2) y; Rr 
SOU >AA>ZZ = Tagh x (1 2 er eal a ar 2 .m? +im,T 
MA MA MA i TIVA Mp; Mht h; 
+ Y (0105/2v) Ro; Rit x (dov,/2v) Rai Ri 
: 4m2, — m. iMh, Lh; i : 4m, — m. impp, 
2 2 2 4 
m2, me, mi,  3mby (02/2) Ry; Ri 
< ov > = 1 1 | x 2 
TUCSWMCS pr m2, ( iy må ) £ 4m?, — m; c impp, 
+ 3 (0105/2v) Ro; Ri j y (dov,/2v) Rai Ri 
; 4m, — m; +imnTn, | | : 4m?, — my. + impp, 
e Top-right 
2 ( 2y )R (R; 929. R; jy + R; d2y 
1 1 må pardon 2i il i2 Vo 139 Vs 
< ov > = 
VV—AA 327 mi, mi, : Am, — my. - impp, 
2 2 2 4 
må mi, m, 3my (02/2) Ry; Ri 
< > = 1 x 
OV 7AA3VV-—- 158a m3 ( m må ) | |: 4m? — m; + impp, 
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(0105/2v) Roi Ria 


(dvs /2v) Rai Ri 


e Bottom-left 
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2 
4m4 


& 


2 . 
— mp, + mn, Vn; 


4m, — m, - imp Tp, 


Í 


2 B 
1 1 m (aive) 
< ov PV shine = dened A y| RA: 
"my my m, — mọ 
2 F 2 2 
p ERG! ma, | (820/2)? Rul, j Ph. | (0919/2). Rai 
AA-hihi ~ Jn m? ni. | me. msn? | "M pu gau 
A A hy A A ha A 
2 
- my. |(dav,/2)? Rai 
| mA | Mhs mA 
e Bottom-right 
2 
1 1 mi. |93 
"E y > deagull = 1 i Rà. 
VV—ohihi 327 mi, mi. 2 2i 
1 1 mi | 62|” mj, [51 vo 
"T j; Seagull = 1 1 R2. | 2 Rà. 
AA-hihi 327 m mA 2 li m2, 2 y 2i 
2 2 
ha do Us 2 
= — 5| m. 
* mà 2v 3i 
B Boltzman equations 
The coupled boltzman equations could be written as [3]: 
dYa 1.32g’ M4Mp 2 2 
: € 
"Pp m ( (ovr) aac (Ya = (1) 
A 
y 
+ (Opel) AASVV (v4- ay n nc)) (B.1) 
(Y 
dYy 1.32917 My Mp 
dry = 3 (OUrel) VV sx € (Y? = po 
V 
yrs 
— (OUrel) AA+VV (v4 ay adj) (B.2) 


for M4 > My. 


(vy 
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Similarly, for My > Ma, one has 
dYy 1.3291 My Mp 
2 


(rove = (mr 


dry V 
(gy 
+ (Gtra)vv-AA (vw — (Yea YA ; (B.3) 
dYa 1.32917 ° M4 M 
dz = E. P (OUrel) AA X X (YÀ = Gay) 
yor 2 
— (0vret)vv Aa (Yv — l v) y? : (B.4) 
(Ya) 
Here, Mp = 2.44 x 10!8 GeV is the reduced Planck mass, and g, is the degrees of freedom 
parameter. Two dimensionless variables Y; j relate with number density through Y; j = Tij ; 
Tij = T [66], and s and T are the entropy density and temperature of the Universe. After 


solving the coupled Eqs. (B.1), (B.2), (B.3), and (B.4), one gets the values of Y4 and Yy 
at present temperature To. 


C One loop £ functions 


The one-loop renormalization group equations (RGEs), which we used to analyze the vac- 
uum stability problem, are given by 


dX 1 


dogm" 16:2/X (C.1) 

with one-loop 6-functions fx, 
Be = 5 (ddp — 99262 + 45? + 38s + 2Amô2) , (C.2) 
Bsa = 5 (dð — 3g265 — 9926» + 129265 + 462 + 1282A + 2ô1 Am) , (C.3) 
NER 5 (1003 +62 462), (C.4) 
Bin = 5 (6162 — 391 Àm — 993Am + 12g? Am — 995Am + 12m + 8AZ, + 3AmAg) ,(C.5) 
Bx = —A(3gi + 995 — 12g? — 24A) + ah + sol + 93) — 6g; 4 S | m , (C6) 
B 5 (908 + 267 + AM, — 1893Ag + 93) , (C) 
Bg, = - S : (C.8) 
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